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摘 要 : 黄土 沟谷 的 发 育 演化 进程 研究 对 于 发 气 沟 谷 形成 与 土壤 侵蚀 的 内 在 机 理 ,指导 黄土 地 区 生态 修复 有 深刻 
的 指导 意义 。 基 于 沟谷 水 平 空间 、 气 候 、 沟 谷 点 线 面 特征 选取 10 个 沟谷 发 育 程度 量化 指标 ,利用 组 合 炉 权 RSR 法 优 


异 的 多 因子 综合 性 分 档 评价 能 力 建立 董 志 震 沟谷 发 育 阶段 分 级 评价 模型 ,探究 黄土 沟谷 的 发 育 演化 及 分 布 规律 。 


人 研究 发 现 :(1) AdETYABBSOGEHTSEH- EUSEB ONDVID JI SAE EY 16.08% ,其 次 为 沟谷 密度 (15.621% ) , 两 者 
是 判断 沟谷 发 育 程度 的 重要 指标 。(2) FETA AIA RSR 法 将 董 志 藉 82 条 沟谷 发 育 程 度 分 为 幼年 期 青年 期 壮年 


期 和 老年 期 ,壮年 期 和 老年 期 沟谷 面积 占 比 为 88.48%， 


E 


ES ZR REUS TERRE A TEEB ,对 于 震 面 侵蚀 强烈 。(3) 


沟谷 发 育 程度 与 其 下 伏 古 地 貌 \ 构 造 地 质 和 黄土 厚度 有 着 高 度 相关 性 。(4) 模型 方差 一 致 性 检验 显著 性 水 平 已 < 


0.001 ,线性 回归 拟 合 优 度 尼 为 0.986 ,评价 模 型 表现 良好 。 


关键 词 : 黄土 ;量化 指标 体系 ; 组 合 炉 权 RSR 法 ; 沟谷 发 育 分 级 评价 ; SOR 


沟谷 作为 黄土 地 区 典型 的 侵 刨 地貌 单 元 ,在 内 
DEJI BER AE DA ZRTER] PARIBUS ei Te. TAS AY 
发 育 程 度 存在 显著 的 差异 性 ,研究 沟谷 发 育 的 发 育 
规律 ,帮助 更 清楚 的 了 解 黄土 沟谷 有 组 织 .系统 的 
发 育 和 演化 ,综合 性 的 量化 指标 也 利于 多 角度 认识 
沟谷 的 纵向 、 横 向 和 水 平方 向 的 特征 ,从 而 帮助 精 
VE EL RAI SK tk E ^] DIS TR Dale > fs 53 PE RJ P 
SENT Az ER PR Hil. AS SCAE TAL MAR RSR 对 董 志 
WRT AY BO EAT aD EC Hf ,帮助 我 们 更 正确 清晰 的 
认识 沟谷 的 发 育 ,为 针对 性 的 减缓 水 土 流失 和 生活 
生产 的 安全 提供 参考 依据 。 

黄土 沟谷 评价 研究 主要 以 侵蚀 类 型 分 区 以 及 
地 貌 类 型 分 类 为 主 "。 基 于 侵蚀 类 型 可 分 为 宏观 的 
黄土 沟谷 侵蚀 分 区 以 及 微观 的 沟谷 侵蚀 发 育 进 
程 研 究 -。 而 从 地 貌 类 型 的 角度 出 发 ,科学 的 沟谷 
地 貌 认 知 及 其 分 类 体系 是 探讨 黄土 沟谷 发 育 过 程 
的 重要 基础 。 长 期 以 来 ,在 DEM 数 据 的 支撑 下 ， 
TUERI" CET JR EJUS "和 景 可 "… 基 
于 长 度 .宽度 .密度 .分支 比 网 络 结构 (如 树枝 状 结 
构 ) 等 不 同 指标 对 沟谷 的 形态 特征 .空间 格局 ,分布 
规律 进行 了 深入 的 研究 ,其 主要 人 研究 了 沟谷 在 水 平 
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方向 上 的 特征 以 及 分 异 情 况 。 何 十 等“\ 陈 传 康 "*、 
严 宝 文 等 近 和 页 爽 等 上 依据 定性 描述 或 半 定 量 的 
方法 对 沟谷 形态 进行 表述 ,从 不 同 的 角度 ,微观 和 
宏观 不 同 尺度 对 黄土 沟谷 体系 进行 了 研究 ,基本 上 
将 沟谷 系统 分 为 细 沟 、 浅 沟 . 切 沟 、 冲 沟 . 干 沟 和 河 
沟 等 类 型 。 选 择 合适 的 沟谷 发 育 量化 评价 参考 指 
标 是 综合 建 模 分 级 评价 的 关键 因素 ,目前 对 于 沟谷 
量化 评价 指标 的 研究 已 经 取得 了 一 定 的 成 果 , 但 大 
都 忽略 了 沟谷 的 层次 性 以 及 系统 性 。 严 宝 文 等 0 
利用 高 程 积分 曲线 模型 将 单一 沟谷 分 为 幼年 期 EE 
年 期 和 老年 期 ; 李 晨 瑞 拉 选取 多 指标 建立 了 有 层 
次 、 系统 性 的 评价 体系 将 黄土 高 原 沟 谷 发 育 程度 分 
为 了 8 个 级 别 ,考虑 了 沟谷 的 系统 性 和 层次 性 来 选 
取 量 化 指标 ,但 未 将 沟谷 集合 作为 评价 对 象 。 熊 礼 
阳 中 认为 流域 由 若干 条 沟谷 组 成 ,各 条 沟谷 的 发 育 
都 不 是 独立 的 ,其 物质 能 量 ,梯度 变化 等 通过 径流 
点 的 串联 来 实现 有 机 连接 ,将 一 个 沟谷 集合 作为 研 
究 对 象 , 但 未 建立 综合 量化 评价 模型 。 

本 文 从 沟谷 点 线 面 特征 、 气 候 和 水 平 空 间 角度 
系统 性 选取 评价 指标 ,并 基于 组 合 炉 权 RSR 法 建立 
综合 评价 模型 ,将 董 志 需 82 个 沟谷 集合 划分 为 4 个 
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^e TELS Be BT Bt , W [125 ER S T FRI JEJE EDU 
考 依据 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

研究 区 位 于 甘肃 省 东部 庆 阳 市 境内 (图 1), 包 
括 33 个 乡镇 ,南北 纵向 跨度 100 km, 东 西 横向 跨度 
46 km, 面积 约 2835 km ,其 中 , Jii TT TALK 968.48 km’, 
Fe TE RS GSC IA] — B HE ,号 称 "“ 天 下 黄土 第 一 
JR". WEIS MRA Bi HEA PA UR , 土壤 为 黑 
3*5 TRUE Zi EE PRA, JEEE E 
为 主 ,厚度 在 50~200 m YD BIE ,边缘 较 薄 。 黑 坊 
土 有 机 质 含 量 丰 富 ,土壤 较为 肥沃 , 黄 绵 土 虽然 有 
机 质 含量 较 低 ,但 可 耕 性 良好 ,是 卉 面 地 势 平 坦 , 人 研 
究 区 土地 资源 优越 ,农业 发 达 。 人 研究 区 呈现 中 间 稍 
高 ,南北 低 的 地 势 格 局 ,腹地 西峰 区 附近 海拔 1400 
m, 北 端 驿 马 镇 一 带 海拔 1300 m ,南部 和 盛 镇 一 带 海 
45.1200 m, 卉 面 周围 受到 冲 沟 切 蚀 呈 树 校 状 和 锯 
TR LIA AE ABE pU TRU , E ZR RC t 
fi. WRK E Ts RK , WANAE 
相间 ,延伸 各 异 ,沟谷 发 育 具有 高 度 可 研 性 ” 。 
1.2 数据 与 方法 
1.2.1 研究 思路 ”本 研究 分 为 3 个 阶段 :(1) 原始 数 
据 的 准备 与 解 译 ; 本 人 研究 的 DEM Landsat 8 谷歌 高 
清 卫 星 影像 (分 状 率 1 m) 、 降 雨量 以 及 路 网 均 为 


2020 年 数据 ;基于 ArcGIS 10.3 软件 ,采用 人 机 交互 
目 视 与 AreToolbox 进行 解 译 。(2) 选取 沟谷 点 线 面 
特征 气候 和 水 平 空间 指标 进行 单一 指标 量化 评 
价 。(3) SEF 4H AP REX RSR JEM GERA AZ E 
行 分 级 评价 ,并 验证 模型 的 准确 性 。 

12.2 参数 选取 ADAH pm, [x 
据 科学 性 原则 .敏感 性 原则 和 可 实现 性 原则 选取 了 
10 个 参考 指标 ,包含 了 正 癌 因子 :边界 维 数 、 沟 谷 密 
E ,地表 切 割 深度 .平均 坡度 和 月 平均 降雨 ,其 量化 
结果 值 越 大 代表 沟谷 发 育 程度 越 高 (如 表 1); 负 问 
因子 : 主 沟 文 沟 比 .土地 利用 综合 指标 、 主 沟 纵 降 
比 、 归 一 化 植被 指数 (NDVI) 和 圆 度 率 , 其 数值 越 小 ， 
代表 发 育 程 度 越 高 ( 表 2)。 这 些 参考 指标 从 气候 、 
水 平 空间 ,沟谷 点 线 面 特征 等 不 同 的 角度 对 沟谷 的 
发 育 程 度 进行 评价 。 

上 述 指标 均 通 过 ArcGIS 10.3 软件 进行 计算 ,得 
SSO EET EARN HA. Wa eT 
表格 显示 分 区 统计 工具 (Zonal Statistics As Table ) Zt 
计 各 量化 指标 数值 在 各 个 沟谷 沟 沿线 区 域内 的 平 
均值 ,以 该 平均 值 作 为 沟谷 发 育 程度 的 评价 依据 。 
123 Maik 考虑 到 不 同 指标 对 于 沟谷 发 育 程 度 
PAF IRAP LE DE EE, RSE RA AO A RE TB on HY AX 
Eo AGA HS SER EAR ET hy EE A) OR f 
定 其 客观 权重 。 其 求解 步骤 如 下 : 

(1) 数据 标准 化 

首先 将 各 个 指标 进行 去 量 纲 化 处 理 。 假 设 给 


图 1 研究 区 地 理 位 置 
Fig. 1 Location of the study area 
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R1 沟谷 发 育 评价 正 向 指标 
Tab.1 Positive indicators for gully development evaluation 
指标 名 称 计算 公式 解释 指标 意义 
边界 维 数 pa 2(In N- C) M HA LER SIN DH AREE SC ”反映 水 平方 向 上 沟谷 空间 展 布 的 复杂 程度 和 沟谷 的 
InM 为 沟谷 周 长 发 育 程 度 。 指 标 值 越 大 ,说 明 局 部 沟谷 与 整体 沟谷 
网 络 的 相似 程度 越 高 ,沟谷 形态 结构 越 复杂 ,发 育 程 
沟谷 密度 5 DIARRE; / 为 流域 内 各 条 沟谷 ”反映 沟 头 前 进 及 沟谷 侵蚀 发 育 程度 的 重要 指标 。 指 
D,= a 的 长 度 ;4 为 流域 面积 标 值 越 大 ,说 明 沟 头 前 进程 度 越 大 , 即 沟谷 生长 延伸 
4 程度 越 大 ,发 育 程度 越 高 。 
地 表 切 市 深度 G= Hua Hus G, 为 地 表 切 制 深 度 ; 及 .为 地 面 某 点 ”直观 反映 了 沟谷 下 切 发 育 程度 和 地 表 被 侵蚀 切割 的 
一 定 邻 域 范围 内 的 平 n 高 程 ; Ha 为， 程度。 指标 值 越 大 ,地 表 侵 蚀 切割 越剧 烈 ,沟谷 下 切 
该 点 同一 邻 域 范围 内 的 最 低 高 程 发 育 程 度 越 高 。 
平均 坡度 slope = arctan [ f +f? x 180. 1go 及 为 zx 方向 的 高 程 变化 率 ; f 为 y 方 向 ”在 一 定 程度 上 反映 沟谷 地 貌 的 发 育 程 度 。 沟 谷 发 育 
的 高 程 变 化 率 初期 ,沟谷 发 育 以 下 切 为 主 ,指标 值 逐渐 增 大 , 随 着 
沟谷 发 育 程度 的 增加 ,沟谷 在 深切 的 同时 也 在 横向 
展 宽 ,指标 值 反而 有 所 下 降 。 
月 平均 降雨 PET R, 为 月 平均 降雨 ; R, 为 年 总 降雨 气候 影响 因素 之 一 ,降雨 的 多 少 一 定 程度 上 影响 水 
12 流 侵蚀 的 剧烈 程度 。 
R2 沟谷 发 育 评价 负 向 指标 
Tab. 2 Negative indicators for gully development evaluation 
指标 名 称 计算 公式 解释 指标 意义 
主 沟 支 沟 比 2b ROB SES XC EU; L 为 主 沟谷 长 度 ; L 从 沟谷 长 度 的 角度 反映 沟谷 系统 中 文 沟 的 发 育 程度 。 指 
L 为 沟谷 总 长 度 标 值 越 小 ,沟谷 系统 中 支 沟 发 育 越 明显 ,说 明 沟 谷 系统 整 
体 发 育 程度 越 高 。 
土地 利用 程度 综 LL, 为 土地 利用 程度 综合 指数 ;4 为 第 i 是 针对 一 个 具体 年 份 内 所 有 土地 利用 类 型 整体 反映 的 土 


L,- 100x Y Ax C, 


合 指标 
主 沟 纵 比 降 nc 
i 
归 一 化 植被 指数 _ NR-R 
EE NDVI WARTR 
pines Nu 
mG 
定 了 p 个 指标 : 


X,={x,, X5. Xa. 777 


假设 各 指标 标准 化 后 的 值 为 Y,Y,, Y, Y, , 


那么 : 
X; 一 min(X,) 
max(X;) 一 min(X;) 


ij 


-" max(X;) = 
i max(X,) - min(X;) 


级 的 土地 利用 程度 分 级 指数 ;C; 为 第 i 地 利用 集约 化 程度 。 通 过 对 土地 利用 程度 的 分 级 ,量化 人 
级 土地 利用 程度 分 级 面积 百分比 类 活动 对 于 土地 系统 的 影响 程度 ,定量 地 描述 该 地 区 土地 
利用 的 综合 水 平和 变化 趋势 ,是 衡量 指定 区 域内 土地 利用 
情况 的 一 个 重要 指标 。 
J 为 主 沟 道 纵 比 降 , H 为 主 沟 道 出 水 口 反映 沟谷 整体 侵蚀 势能 ,是 描述 沟谷 在 垂直 方向 上 侵蚀 发 
点 的 高 程 ; Hy 为 主 沟 道 沟 头 点 的 高 程 ; 育 程度 的 重要 指标 。 指 标 值 越 小 ,说 明 沟谷 的 侵蚀 能 量 越 
1 为 主 沟 道 长 度 乓 ,沟谷 越 趋 近 于 稳定 的 宽 谷 ,发 育 程度 越 高 。 
NIR 为 近 红 外 波段 的 反射 值 ;R 为 红 光 反映 土地 覆盖 植被 状况 的 一 种 遥感 指标 ,定义 为 近 红 外 通 
波段 的 反射 值 道 与 可 见 光 通道 反射 率 之 差 与 之 和 的 商 。 
e WARR; S, 为 沟谷 面积 ; C, 为 沟 沿 反应 水 平方 向 上 沟谷 接近 圆 的 程度 ,其 值 越 接 近 1, 则 越 
线 长 度 接近 圆 ,其 发 育 程度 越 低 。 
P= : , (i=l, ,1;]=1, ^: (4) 
， xj (1) Y. 
ij 


CIEIBHRER ) 


( 负 向 指标 ) 


(2) 求 各 指标 在 各 方案 下 的 比值 


第 7 项 指标 在 第 
计算 该 指标 的 变异 大 小 。 


p 


(2) 


(3) 


第 i 个 方案 中 占 该 指标 的 比重 , 即 


(3) ae TES rt n Ad 


A d fei e YE n fe TRE DAI a. X ,一 组 数据 的 信息 
WIN: 
Ez -In(n) S P nP, (5) 
式 中 : ERO. P,=0 dux 已 =0。 
(4) — 
w,= (j= 1,2, ---,p) (6) 


YE 
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式 中 ;为 指标 个 数 , 即 Ep. 
1.24 秩 和 上 比 综合 评价 法 秩 和 比 综合 评价 法 
(RSR 法 ) 是 一 种 将 古典 参数 统计 和 近代 非 参 数 统 
计 进 行 结合 ,并 融 其 各 自 优 点 于 一 身 的 统计 分 析 方 
法 。 秩 和 比 综合 法 根据 计算 的 无 量 纲 统计 量 RSR 
进行 分 档 评 价 ,避免 一 般 评价 方法 局 限 性 ,运用 多 
个 指标 对 多 个 单位 进行 评价 ,具有 综合 性 和 系统 
性 ,使 评价 变 得 客观 准确 。 

其 计算 步骤 如 下 : 

(1) 列 出 原始 数据 表 

列 出 原始 数据 表 。 假 设 有 nn 个 待 评价 样本 ,p 项 
评价 指标 ,形成 原始 指标 数据 矩阵 : 
l (7) 
AP: w 表示 第 ;个 样本 第 1 项 指标 的 数值 。 

(2) 计算 秩 值 

对 评价 指标 按 其 指标 值 的 大 小 进行 排序 ,得 到 
秩 次 RR, 用 秩 次 RR 来 代替 原来 的 评价 指标 值 。 人 研究 
选用 非 整 秩 法 进行 编 秩 ,结果 建立 各 指标 的 秩 次 数 


JEPE. 
Exc (8) 
R, Pr R, 


式 中 : Ry 表示 第 :个 样本 第 7 项 评价 指标 的 秩 次 。 
对 于 正 向 指标 : 


X- 


=X. 


X; min 
R; 7 : "o 7 U X max =Å ain 9 
对 于 负 向 指标 : 
R =1 +- oma! (10) 
i 7 an i Pm - A ois 
RP: X. =max(X,, Xj «Xj; X = min(X,. X, 
a A 
(3) 计算 秩 和 比 RSR 值 及 排名 


在 一 个 n 行 p 列 矩阵 中 ,RSR 的 计算 公式 为 : 
"EE. 
RSR, = Ta 
AP: i=1,…,n; j=l, p, Ry 表示 第 i 行 第 j 列 
元 素 的 秩 。 


赋予 不 同 指标 权重 则 计算 加 权 秩 和 比 为 
WRSR ,其 公式 为 : 


(11) 


(12) 


式 中 : W, 表示 第 个 评价 指标 的 权重 ,满足 SW, = 1 。 


(4) 确定 RSR 的 分 布 

将 RSR 值 按 从 小 到 大 顺序 排列 ,计算 出 各 组 频 
数 、 累 计 频 数 、 秩 次 R 以 及 平均 秩 次 户 以 及 向 下 累 
计 频 率 R nx 100% ,最 后 一 项 用 (1-1/4n)x100% 修 
正 ,查询 百分数 与 概率 单位 对 照 表 得 到 概率 单位 


Probit 值 。 
(5) 基于 SPSS 软件 计算 回归 方程 ,得 到 拟 合 
方程 。 


(6) 据 各 分 档 排 序 情况 下 概率 单位 Probit 值 , 按 
照 最 佳 分 档 原则 对 评价 对 象 进行 分 档 归 类 。 分 档 
数 由 研究 者 根据 实际 情况 决定 。 
12.5 模型 准确 性 验证 方差 齐 性 检验 (Homoge- 
neity of Variance Test) 是 数理 统计 学 中 检查 不 同样 
本 的 总 体 方差 是 否 相 同 的 一 种 方法 。 其 基本 原理 
是 先 对 总 体 的 特征 作出 某 种 假设 ,然后 通过 抽样 研 
究 的 统计 推理 ,对 此 假设 应 该 被 拒绝 还 是 接受 作出 
HEMT. TFA TIAA Hartley 检验、Bartlett 检 验 、 修 正 
的 Bartlett 检 验 。 研 究 进行 方差 一 致 性 检验 (Bartlett 
检验 ) 检 验 模型 的 准确 性 。 

线性 回归 分 析 宇 所 得 的 尼 代 表 曲 线 的 拟 合 程 
E ,其 值 越 接近 于 1 代表 拟 合 程度 越 好 。 方 差 一 至 
性 检验 (Bartlett 检验 ) 要 求 各 档 差 异 具 有 统计 学 意 
Xo VIF 值 用 于 验证 模型 是 否 存在 共 线 性 问题 , VIF 
值 小 于 10 或 5, 则 模型 各 指标 之 间 不 存在 共 线 性 


问题 。 
2 结果 与 分 析 


2.1 数据 解 译 与 指标 计算 

本 文 使 用 AreGIS 10.3 软件 基于 2020 年 的 谷歌 
高 清 卫 星 影 像 ( 分 辨 率 1 m) 目 视 解 译 董 志 抓 地 区 82 
条 沟谷 水 系 网 络 以 及 沟谷 沟 沿 线 (图 2)。 

黄土 高 原 上 沟谷 存在 纵向 .横向 和 垂 向 3 种 侵 
蚀 方 式 , 沟 头 的 溯源 侵蚀 , 沟 坡 的 侧 蚀 以 及 下 切 侵 
蚀 。 选 取 的 10 个 指标 中 ,地 表 切 割 深度 以 及 主 沟 纵 
降 比 反映 了 整个 沟谷 垂 向 的 下 切 侵蚀 的 程度 ; 圆 度 
率 .边界 维 数 在 一 定 程度 上 反映 了 沟谷 水 平方 向 上 
的 发 育 水 平 。 其 他 指标 从 气候 、 水 平方 向 的 几何 复 
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3 作为 黄土 地 区 主要 水 源 , 决 定 了 径流 量 的 大 小 , 整 

Tics Sk KETI 2E EIN, ,沟谷 平均 月 降雨 量 均 在 
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图 2 董 志 震 沟谷 原始 数据 解 译 
Fig.2 Interpretation of raw data of Dongzhi Loess Plateau Gully 


杂 度 和 植被 等 角度 进行 评价 。 

使 用 ArcGIS 10.3 对 所 选取 的 指标 计算 得 到 各 
沟谷 各 指标 的 平均 值 作为 沟谷 发 育 程 度 评判 依 
据 。 将 各 指标 数值 赋值 于 各 沟谷 ,得 到 沟谷 发 育 程 
度 分 异 图 (图 3, 图 4)。 正 向 指标 中 ,平均 坡度 在 一 
定 程度 反映 沟谷 下 切 的 程度 ,其 值 在 东部 及 北部 最 
大 达到 24.5 ,西部 南部 则 为 18 左右 ;边界 维 数 从 水 
平 空 间 描 述 沟谷 的 几何 形状 的 复杂 度 来 评价 沟谷 
的 发 育 程 度 , 北 部 和 东部 地 区 边界 维 数 在 1.3~1.35， 
西部 和 南部 在 1.15$~1.25; 地 表 切 制 深度 反应 沟谷 下 
切 发 育 程度 和 地 表 被 侵蚀 切割 的 程度 ,东部 地 区 切 


(b) 


>z 


地 表 切 割 


平均 坡度 边界 维 数 深度 /m 


uH 17.5~19 E 1.15-12 wm 90-130 
© 19-20.5 = 1.2-1.25 = 130-170 
ES 20.5-23 ES 1.25-1.3 B® 170-210 
EU 23-24.5 ES 1.3~1.35 um 210-260 


0 10km 0 10km 0 10km 


45.5-49.5 mm; 沟谷 密度 反映 沟 头 前 进 及 沟谷 侵蚀 
发 育 程度 ,东部 和 西部 值 较 大 ,集中 在 5~6 km* km", 
北部 集中 在 3~5 km*km*。 负 向 指标 中 , 主 沟 支 沟 比 
从 沟谷 长 度 的 角度 反映 沟谷 系统 中 支 沟 的 发 育 程 
度 , 东 部 和 北部 值 较 小 为 0~0.1, 西 部 和 南部 在 0.1~ 
0.5; 主 沟 纵 降 比 在 垂直 方向 上 反映 侵蚀 发 育 程度 ， 
东部 和 北部 大 部 分 地 区 处 于 0~0.2 ,西南 部 和 南部 处 
于 0.2~0.8; 圆 度 率 反映 水 平方 向 上 接近 圆 的 程度 ， 
东部 和 西北 部 沟谷 沟 沿 线 最 不 接近 圆 ,在 0~0.12, 北 
部 在 0.24~0.48 最 接近 圆 ;NDVI 反 映 植被 覆盖 的 状 
况 ,其 值 越 高 代表 植被 覆盖 越 好 ,东部 集中 于 0.28~ 
0.33 ,北部 集中 于 0.23~0.28 ,南部 和 西部 部 分 地 区 在 
0.33-0.4; 土地 利用 程度 综合 指标 评价 了 土地 利用 
的 集约 化 程度 ,也 在 一 定 程度 上 反映 了 植被 的 分 布 
趋势 ,东部 和 西部 绝 大 部 分 都 在 240~275 ,北部 则 落 
在 205~240。 

正 向 指标 东部 和 北部 计算 值 比 西部 和 南部 较 
大 , 负 向 指标 显示 东部 和 北部 计算 值 比 西部 和 南部 
BUR, GS ERAI TE in, IISA 2 Fe BOE ACA 
致 为 偏向 东 侧 及 北 侧 沟 谷 发 育 程度 较 高 ,偏向 西 侧 
及 南 侧 沟 谷 发 育 程 度 较 低 。 

22 董 志 震 沟谷 发 育 程度 分 级 评价 

本 研究 中 ,使 用 炉 权 法 计算 出 所 选取 的 10 个 指 
标的 权重 (图 5)。 其 中 , NDVI 所 占 权 重 最 大 为 
16.08% ,平均 坡度 所 占 权 重 最 小 为 4.823% ,其 余 指 
标 权 重 分 别 为 分 形 维 数 的 权重 为 13.222% 、 沟 谷 密 
度 的 权重 为 15.621% 地表 切割 深度 的 权重 为 


(d) 


À 


(e) 


À 


月 平均 沟谷 密度 
降雨 量 / Kkm*km?) 


EB 45-46 EB 3-4 


E 46-47 m 4-5 
Ed 47-48 ES 5-6 
轩 48-49 6~7 
0 10km 0 10km 


图 3 沟谷 发 育 正 向 指标 赋值 


Fig. 3 Gully development positive indicator assignment 
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(a) 


m 


土地 利用 Md 
NDVI 综合 指数 
ES 0.18~0.23 ES 170-205 
UC 0.23~0.28 © 205-240 
BW 0.28~0.33  240~275 
0.33~0.4 275~310 
0 10km 0 10km 


图 4 沟谷 发 育 负 向 指标 赋值 


Fig.4 Gully development negative indicator assignment 


主 沟 支 
沟 比 圆 度 率 
= 0~0.1 EB 0-0.12 
= 0.1-0.2 = 0.12-0.24 
ES 0.2-0.3 ES 0.24-0.36 
0.3-0.4 0.6-0.8 EB 0.36-0.48 
0 10km 0 10km 0 10km 
— 边界 维 数 
| AV J 
综合 指数 eae! 
月 平均 
NDVI 了 南明 
Xxx PIE 
a gun 
圆 度 率 平均 坡度 
主 沟 纵 降 比 


图 5 沟谷 发 育 程度 评价 指标 权重 
Fig.5 Weighting of indicators for evaluating the degree of 
gully development 


8.997% 、 降 雨 的 权重 为 6.372% 、 圆 度 率 的 权重 为 
7.766% 、 主 沟 纵 比 降 (11000) 的 权重 为 8.835% 、 主 
沟 支 沟 比 的 权重 为 7.609% .土地 利用 综合 指数 的 权 
重 为 10.676%。 

计算 得 到 的 每 条 沟谷 各 指标 的 RSR 值 ,利用 插 
值 法 参考 百分数 与 概率 单位 对 照 表 ,得 到 其 相对 应 
的 Probit 值 。 其 中 RSR 临界 值 ( 拟 合 值 ) 为 根据 对 应 
的 Probit 值 代入 拟 合 方程 中 求解 得 到 。 依 据 每 个 沟 
谷 的 拟 合 RSR 值 对 应 的 Probit 值 从 大 到 小 进行 排列 


分 档 评价 , 百 分 位 数 临 界 值 和 Probit 临界 值 根 据 分 
档 水 平 数量 而 变化 ,该 两 项 是 固定 值 且 完全 一 一 对 
应 ( 表 3)。 

基于 组 合 箭 权 RSR ZEE A AT 4 
个 阶段 :幼年 期 青年 期 壮年 期 和 老年 期 (图 6, 表 
3)。 沟 谷 的 溯源 侵蚀 主要 将 沟谷 拉 长 ,下 刨 则 是 加 
深 沟谷 , 侧 蚀 从 水 平方 向 上 展 宽 沟 谷 。 董 志 基 北 部 
沟谷 多 为 壮年 期 ,面积 占 比 达 49.06% ,水 平方 向 上 
的 发 育 空间 十 分 受 限 , 侧 蚀 及 溯源 侵蚀 作用 较 弱 ， 
但 通过 地 表 切 割 深度 和 主 沟 支 沟 比 判断 ,整体 在 垂 
向 上 的 下 蚀 作 用 仍 有 较 大 的 发 育 空间 。 幼 年 期 沟 
谷 占 比 0.63% ,处 于 发 育 的 初期 ,各 项 指标 都 有 很 大 
的 发 育 空间 。 青 年 期 占 比 10.89% ,沟谷 在 水 平方 向 
和 垂直 方向 上 均 有 较 大 的 发 育 空间 和 发 育 潜力 。 老 
年 期 占 比 39.42% ,沟谷 各 项 指标 都 表明 沟谷 水 平 扩 
张 和 下 切 发 育 较 完备 ,整体 沟谷 形态 较为 复杂 *””。 
2.3 线性 回归 分 析 及 方差 一 致 性 检验 

基于 SPSS 软件 进行 线性 回归 分 析 和 方差 一 致 
性 检验 ,模型 表现 良好 ( 表 4):(1) 对 于 变量 共 线 性 
表现 ,VIF 全 部 小 于 10, 因 此 ,模型 没有 多 重 共 线 性 
问题 ;(2) 线性 依赖 于 其 未 知 参数 的 模型 比 非 线 性 
依赖 于 其 未 知 参数 的 模型 更 容易 拟 合 ,而 且 产 生 估 


H 


A3 董 志 震 沟谷 发 育 阶段 分 级 
Tab.3 Stages of development of the Dongzhi Loess Plateau Gully classification 


发 育 阶段 百 分 位 临界 值 Probit 秩 和 比 临 界 值 ( 拟 合 值 ) 频数 apu 面积 占 比 
幼年 期 <6.681 <3.5 <0.3106 5 10.04 0.63% 
青年 期 6.681~50 3.5~5 0.3106~0.4692 36 177.61 10.89% 
壮年 期 50~93.319 5-6.5 0.4692~0.6277 35 799.58 49.06% 
老年 期 > 93.319 > 6.5 > 0.6277 6 642.44 39.42% 
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图 6 董 志 震 沟谷 发 育 程度 分 异 
Fig.6 Variation in the degree of development of the Dongzhi 
Loess Plateau Gully 


计 的 统计 特性 也 更 容易 确定 。 沟 谷 发 育 评价 指标 
的 RSR 值 拟 合 方程 y=-0.059+0.106xProbit ,模型 的 
拟 合 优 度 尼 为 0.986, 模 型 表现 较为 优秀 ; (3) 从 方 
差 一 致 性 检验 (Bartlett 检 验 ) 的 结果 分 析 可 以 得 到 ， 
显著 性 水 平 P<0.001, 水 平 呈 现 显著 性 ,拒绝 了 回 
归 系 数 为 0 的 原 假设 ,模型 基本 满足 要 求 。 


3 讨论 


基于 烂 权 法 计算 得 到 权重 最 高 的 两 个 指标 为 
NDVI(16.0896 ) 和 沟谷 密度 (15.621% ) ,其 大 小 是 判 
斯 沟谷 发 育 程度 的 重要 依据 。 植 被 极 大 程度 上 影 
响 了 土壤 侵蚀 ,其 减少 水 土 流失 的 机 理 主要 为 : (1) 
冠 层 减弱 雨滴 的 一 定 势能 和 动能 ;(2) 落叶 调节 表 
面 径 流 ;(3) 影响 土壤 降水 入 渗 “;(4) 根系 可 提高 
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土 体 抗 侵蚀 和 抗 剪 强度 ,黄土 坡 面 主要 以 均 质 士 固 
结 为 特点 ,植物 根系 与 土壤 间 的 静摩擦 阻力 是 决定 
根系 提高 土壤 抗 剪 强度 增 量 的 主要 因素 ”2 APR 
曙 振 的 研究 表明 随 着 林 龄 的 增长 ,其 根系 对 于 土 体 
抗 剪 强度 提高 越 明显 ( 表 5)。 植 被 减少 水 土 流 失 的 
同时 ,土壤 侵蚀 也 影响 植被 的 发 育 ””。 土 壤 侵 蚀 是 
旱地 和 生态 系统 退化 的 最 重要 驱动 力 ”™, 并 严重 干 
扰 植 物 群 落 发 展 和 植被 演 蔡 过 程 , 从 种 子 形成 开 
始 , 影 响 整 个 生长 阶段 ,对 于 种 子 活力 扩散、 发 芽 、 
幼苗 建立 ,植物 群落 结构 和 植被 空间 分 布 均 造 成 负 
影响 。 例 如 ,土壤 侵蚀 通过 减少 土壤 养分 和 水 分 
来 影响 种 子 的 正常 发 育 ,导致 夭 胎 形 成 失败 和 种 子 
活力 丧失 ,水土 流失 严重 地 区 的 植物 种 子 活 力 明 显 
低 于 土壤 稳定 区 ”…。 植 被 覆盖 于 土壤 侵蚀 之 间 存 
在 互 馈 作用 ,Zhou 等 吕 基 于 气候 变化 .生态 压力 
素 和 人 类 活动 开发 了 一 个 估算 植被 覆盖 和 土壤 侵 
蚀 的 模型 ,土壤 侵蚀 与 植被 覆盖 的 回归 表明 ,一 年 
内 植被 覆盖 度 增 加 1% ,可 以 减少 土壤 侵蚀 456 te 
km .a …, 当 植被 覆盖 率 约 为 32% 时 ,土壤 侵蚀 趋 于 
m. 
EH e 55 0^ i E reg Js DC b ER AA 
度 线性 拟 合 ,得 到 沟谷 密度 与 输 沙 模 数 之 间 存 在 较 
强 的 相关 性 (图 7) ,沟谷 密度 是 反应 沟 头 前 进 及 沟 
谷 侵蚀 发 育 程度 的 重要 指标 ,指标 值 武大 ,说明 沟 
头 前 进程 度 越 大 , 即 沟谷 生长 延伸 程度 越 大 ,发 育 
程度 越 高 。 总 侵蚀 和 非 表 层 土 侵 蚀 主要 受 沟谷 密 
度 变 化 控制 ,沟谷 密度 反 过 来 又 受 地 形 因 子 和 ND- 
VI 的 控制 后 。 沟 谷 密度 的 大 小 受 多 方面 因素 影响 : 
(1) 黄土 高 原 的 地 形 发 育 特 征 受 晚 新 生 代 以 来 的 新 
构造 运动 控制 形成 。 黄 土 高 原 基 本 上 处 于 以 垂直 
方向 差异 性 运动 为 主 的 新 构造 运动 区 域内 ,控制 着 
全 高 原 区 的 地 貌 骨 架 发 育 过 程 和 发 展 方向 ,以 及 水 
系 的 发 育 ;(2) 董 志 藉 地 区 降雨 是 径流 的 主要 来 源 ， 
年 降雨 量 虽然 小 但 集中 于 7 一 9 月 ,降雨 作为 黄土 地 
区 侵蚀 重要 外 力 的 同时 又 促进 了 植被 的 发 育 减 少 


表 4 方差 一 致 性 检验 及 线性 回归 分 析 


Tab. 4 Consistency of variance test and linear regression analysis 


非 标准 化 系数 标准 化 系数 t P VIF R 调整 尼 F 
B 标准 误 Beta 
常数 -0.059 0.007 - -8.173 0.001 
0.986 0.986 5567.72 
Probit 0.106 0.001 0.993 74.617 0.001 1 


注 :8 表示 自 变 量 与 因 变 量 负 相 关 ; Beta 为 标准 化 系数 ;VE 表示 检验 模型 是 否 呈 现 共 线 性 ; 尼 代 表 曲 线 回 归 的 拟 合 程度 ;已 代表 显著 性 水 平 ; 


1 为 7 检验 结果 ;- 表 示 无 数据 。 
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表 5 林 龄 与 根系 提高 土 的 抗 剪 强度 的 关系 
Tab.5 Relationship between the age of the forest and the 
ability of the root system to increase the shear strength of 
the soil 
林 龄 /a 6~10 10~19 20-29 30~40 >40 
提高 土 的 抗 剪 强度 /kPa 0.04 0.08 — 0.13 0.15 0.18 


输 沙 模 数 /10; km? 
a 


Y=728.8x+349.9 
R?—0.63 


1 1 1 1 1 L J 
2 4 6 8 10 12 14 


ARR BE (km: km?) 


图 7 沟谷 密度 与 输 沙 模 数 关系 
Fig. 7 Relationship between gully density and 


sand transport modulus 


水 土 流失 ”。 

董 志 震 沟谷 发 育 程度 各 异 ,其 下 伏 构 造 地 质 、 
黄土 厚度 和 节理 优势 走向 从 根本 上 决定 了 其 差异 
E. ERRATEA RKE PER SR der aa m 
西南 部 , 属 华北 克拉 通 的 次 级 构造 单元 。 地 质 构造 
导致 了 研究 区 在 接受 黄土 堆积 的 同时 构造 侵蚀 也 
在 不 断 加 强 ,其 对 沟谷 侵蚀 发 育 .土壤 侵蚀 等 的 影 
I] PE: C1) 构造 运动 造成 地 面 整体 抬升 与 下 降 
导致 的 侵蚀 基准 面 差异 影响 了 研究 区 侵蚀 空间 条 
件 发 展 ,研究 区 处 于 鄂尔多斯 盆地 的 中 央 古 隆起 南 
部 与 西部 的 边缘 ,北部 为 伊 盟 隆起 ,西部 为 天 环 向 
斜 , 隆 起 构造 内 部 上 升 差异 很 小 ,但 通常 核 部 上 升 
幅度 最 大 ,边缘 上 升 相对 较 小 ,造就 了 研究 区 东 高 
西 低 , 北 高 南 低 的 下 伏 古 地 貌 格局 ,北部 及 东部 侵 
蚀 基 准 面 高 ,沟谷 发 育 程度 也 越 高 (图 8);(2) 研究 
p Pet HA A xL T EE DX E ERE ui A TEN 
厚 , 南 厚 北 薄 的 布局 分 异 ”( 图 9)。 黄 土 堆积 越 厚 ， 
HIRI H R, ARRA TEER” ,同时 ,沟谷 
的 发 育 也 降低 了 沟谷 周围 的 土壤 质量 ,进一步 加 剧 
侵蚀 退化 ” ; (3) 下 伏 基 岩 和 黄土 体内 部 由 于 构造 
地 质 作用 发 育 大 量 的 构造 节理 和 裂 际 。 人 研究 区 存 


CX) Sim 
研究 区 


图 8 董 志 震 沟谷 地 质 构造 
Fig.8 Geological formations of the Dongzhi 


Loess Plateau Gully 


黄土 厚度 /m 
= 50-125 
mm 125-170 
== 170~200 
2200 


0 5km 


图 9 ERRARE EEE 
Fig.9 Loess thickness of the Dongzhi 
Loess Plateau Gully 


202303.10339v1 


chinaXiv 


3 期 1E 


在 3 组 近 正 交 的 节理 系统 :N-S 5; E-W, ENE-WSW 
与 NNW-SSE、NNE-SSW 55 WNW-ESE, ,研究 区 沟谷 
水 系 的 优势 走向 为 N-S 与 E-W NE-SW 与 NW-SE 两 
组 近乎 正 交 的 沟谷 水 系 走向 ,这 与 研究 区 N-S 5 E- 
W All NNE-SSW 5j WNW-ESE 两 组 节理 走向 几乎 一 
致 ,表明 构造 运动 引起 的 节理 裂隙 控制 着 沟谷 水 系 
BAS, 


4 结论 


本 文 基于 遥感 影像 解 译 并 选取 多 角度 描述 沟 
谷 发 育 程 度 的 量化 指标 , ES HAA RSR 模型 对 
董 志 震 沟谷 发 育 程 度 进 行 分 级 评价 ,采用 方差 一 至 
性 检验 和 线性 回归 分 析 验 证 模型 的 可 靠 性 。 最 后 
分 析 董 志 震 沟谷 发 育 的 差异 ,为 开展 针对 性 的 水 土 
保持 提供 参考 依据 。 

(1) 基于 焙 权 法 计算 各 指标 权重 , 归 一 化 植被 
指数 (NDVI) 和 沟谷 密度 是 评价 董 志 韦 沟谷 发 育 程 
度 最 敏感 的 两 个 指标 ,植被 指标 是 沟谷 侵蚀 程度 的 
重要 评判 依据 ,也 是 影响 水 土 流失 的 主要 因素 。 

(2) ETAR RSR 法 建立 沟谷 分 级 评价 模 
型 ,发现 董 志 坟 北部 及 东部 沟谷 的 发 育 程 度 较 高 ， 
对 扰 面 的 蛋 食 严重 ;南部 及 西部 沟谷 密集 ,但 发 育 
程度 不 高 ;处 于 壮年 期 和 老年 期 的 沟谷 面积 占 比 达 
88.48% , 15] 4 Hi ROG] F 3866 EE E GS t p JE F3 R 
RU BEAT ESTA D PR TEJELBE s 

(3) 董 志 卉 沟谷 发 育 的 差异 与 其 构造 地 质 \ 优 
势 节 理 走向 和 黄土 厚度 密 不 可 分 。 董 志 卉 东 高 西 
低 , 北 高 南 低 的 下 伏 古 地 貌 格局 ,北部 及 东部 侵蚀 
基准 面 高 ,黄土 堆积 厚度 也 越 厚 ,沟谷 发 育 程度 越 
高 ;沟谷 水 系 优势 走向 与 N-S 与 E-W 和 NNE-SSW 与 
WNW-ESE 两 组 构造 节理 几乎 一 致 ,表明 节理 裂 际 
控制 着 沟谷 发 育 。 
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Evaluation of Dongzhi Loess Plateau Gully development based on combined 
entropy weight Rank-Sum Ratio method 
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Abstract: As a typical erosion geomorphic unit in the loess areas, gullies continuously erode loess surface under 
the action of internal and external forces, rainfall, and other factors. There are significant differences in the degree 
of development of gullies and valleys. Evaluating the development pattern of gullies and valleys can enhance the 
understanding of their organized, systematic formation and evolution. In addition, comprehensive quantitative 
indexes can facilitate the understanding of the longitudinal, lateral, and horizontal characteristics of gullies and 
valleys from multiple perspectives, which can help in the accurate and effective implementation of gully 
consolidation and loess protection projects, and reduce the constraints on development caused by this fragile 
geological environment. Evaluating the development of loess valleys is also crucial for exploring the intrinsic 
mechanism of valley formation, soil erosion, and guiding the ecological restoration of loess areas. In this study, 
ten quantitative indicators of gully development were selected based on the horizontal spatial, climatic, and point- 
line characteristics of the gullies, and used to establish an indexed evaluation model of the developmental stages 
of Dongzhi Loess Plateau gullies based on the efficient multi-factor comprehensive indexing evaluation capacity 
of the combined entropy- weight rank-sum ratio (RSR) method. Results showed that: (1) The highest weight of 
16.08% could be attributed to the normalized vegetation index (NDVI) based on the entropy weighting method, 
followed by gully density (15.621956), both of which are important indicators for predicting the degree of gully 
development. (2) The combined entropy weighting RSR method could classify the 82 gullies of Dongzhi Loess 
Plateau into juvenile, youth, prime, and old age categories, with the percentage of gully area under the prime and 
old age categories accounting for 88.48% and 88.4896, respectively, while the eastern and northern parts of 
Dongzhi Loess Plateau exhibited more developed features and strong erosion on the loess surface. (3) The degree 
of gully development was highly correlated with the underlying palaeomorphology, tectonic geology, and loess 
thickness. (4) A significant model variance consistency test at P<0.001, with a linear regression R’ fit of 0.986 
was observed, which indicated a good performance of the model. 
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